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Abstract. We present a model of the vertical stratifica- 
tion and the spectra of an irradiated medium under the 
assumption of constant pressure. Such a solution has prop- 
erties intermediate between constant density models and 
hydrostatic equilibrium models, and it may represent a 
flattened configuration of gas clumps accreting onto the 
central black hole. Such a medium develops a hot skin, 
thicker than hydrostatic models, but thinner than con- 
stant density models, under comparable irradiation. The 
range of theoretical values of the index is comparable 
to those from hydrostatic models and both are close to 
the observed values for Seyfert galaxies but lower than in 
quasars. The amount of X-ray Compton reflection is con- 
sistent with the observed range. The characteristic prop- 
erty of the model is a frequently multicomponent iron 
Kq line. 

Key words: Radiative transfer, Accretion disks, Galax- 
ies:active, Galaxies:Seyfert, X-ray:stars, X-rays: galaxies 



1. Introduction 

The intense production of hard X-ray emission is a basic 
property of all active galactic nuclei (AGN), although the 
fraction of this emission in comparison to the total bolo- 
metric luminosity is different in various types of AGN. A 
fraction of this radiation illuminates the relatively cold 
material present there and provides a reflection compo- 
nent apparent in X-ray spectra (Pounds et al. 1990; for a 
review see Mushotzky et al. 1993). The component con- 
sists of an iron Kq, line, usually at 6.4 keV appropriate 
for only weakly ionized iron, an accompanying edge, and 
a Compton "hump" peaking at about 30 keV. 

Above 1 keV the spectral index is equal to 0.7 on av- 
erage, and in a majority of objects the continuum below 1 
keV displays an excess when compared to the extrapola- 
tion of the continuum in the 1-10 keV range (the so-called 
"soft X-ray excess" ) . Recent observations performed with 
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Chandra show that the strength of the excess has proba- 
bly been overestimated in the past. This soft X-ray excess 
may, or may not, be related to the X-ray reprocessing (e.g. 
Gzerny & Zycki 1994, Leighly 1999). 

Above 1 keV the spectrum can be decomposed into 
an underlying power law (with a typical spectral index 
F ~ 1.9) and a reflection component. The strength of the 
reflection component correlates strongly with the slope of 
the power law, and this effect is universal for all accret- 
ing black holes, i.e. AGN and galactic black holes (Ueda 
et al. 1994, Magdziarz et al. 1998, Done & Zycki 1999, 
Zdziarski et al. 1999, Gilfanov et al. 1999, Revnivtsev et 
al. 2001). Frequency-resolved spectroscopy puts additional 
constraints on the formation of the reflection component 
(Gilfanov et al. 2000, Revnivtsev et al. 1999, Zycki 2002) 
but this powerful technique was not applied so far to AGN 
observations. 

The iron line displays an asymmetrical profile with 
an extended red wing and a much smaller blue wing (Nan- 
dra et al. 1997). The presence of the broad red wing in 
the Seyfert composite spectrum has been recently ques- 
tioned (Lubihski & Zdziarski, 2001). However, in the two 
best studied objects, MCG-6-30-15 and NGC 3516, the 
line wing extends down to 4 keV (see Tanaka et al. 1995, 
Iwasawa et al. 1996, Lubihski & Zdziarski 2001 for MCG- 
6-30-15, and Nandra et al. 1999 for NGC 3516). Recent 
high resolution Chandra data show a neutral very narrow 
iron line in most objects (e.g. Yaqoob et al. 2001a; Sam- 
bruna et al. 2001; Kaspi et al. 2001). These results most 
probably mean that the line is multicomponent, as actu- 
ally shown by Yaqoob et al. (2001b). The equivalent width 
of the line seems also to correlate with the spectral slope 
of the power law (Lubihski & Zdziarski, 2001). 

The iron line and the reflection hump have been now 
extensively used for the diagnosis of the geometry of the 
accretion fiow pattern close to the black hole (e.g. Fabian 
et al. 2000, Yu & Lu 2000, Ruszkowski et al. 2000, Abras- 
sart & Dumont 1998, 2000, Karas et al. 2000, Zycki & 
Rozahska 2001, and references therein). 
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Recent intensive monitoring campaigns brought, how- 
ever, more confusion than enUghtening to the emerging 
picture. According to predictions based on a simple anal- 
ysis, a variable X-ray emission should give in response a 
significant variability in the optical/UV band due to the 
absorption of a fraction of incident X-rays. Also both the 
reflected X-ray continuum and the iron line should vary in 
concert, though with a small delay (e.g. Rokaki et al. 1993, 
Clavel et al. 1992, Berkley, Kazanas & Ozik 2000, Kazanas 
& Nayakshin 2001). Recent results of careful monitoring, 
however, showed that the level of the optical/UV response 
is surprisingly different in various sources, the correlation 
between the optical/UV and X-rays is far from perfect or 
is actually absent, the iron line total flux responds only 
weakly or not at all to variable X-rays, and the level of 
reflection in continuum is not positively correlated with 
the strength of the line. 

For instance, NGC 7469 shows variations of the UV 
and X-ray fluxes with similar large amplitudes and a time 
delay of the order of 2 days, the UV leading the X-rays 
in the flux peaks, but the X-ray flux displays in addition 
short term variations not seen in UV (Nandra ct al. 1998). 
The UV flux correlates with T (Nandra et al, 2000). For 
NGC 3516 there is no clear correlation between UV and 
X (Edclson ct al. 2000). For NGC 4051 no optical vari- 
ations were observed when strong X-ray variations were 
seen (Done et al. 1990). The UV flux variations drive the 
optical flux variations with a time delay < 0.2d for NGC 
4151 (Edelson et al. 1996), < 0.15d for NGC 3516 (Edel- 
son et al. 2000), and for NGC 7469 a delay of the order a 
day which increases with the wavelength (Wanders et al. 
1997, Collier et al. 1998). Soft X-ray variations are gen- 
erally larger than hard X-ray ones, but this is not always 
the case (Nandra ct al. 1997). For NGC 3516, there were 
strongly correlated variations with no measurable lag (< 
0.15d, Edelson et al. 2000). 

The ccntroid energy, the intensity, and the equivalent 
width of the Fe line, are variable on short time scales 
(down to ksec). The correlations between the continuum 
and the Fe line have been intensively looked for in a few 
objects with ASCA and RXTE, and now with Chandra, 
revealing a complex behaviour of the line versus the un- 
derlying continuum. Iwasawa et al. (1996) have observed 
rapid variations of the Iron line proflle and intensity in 
1 to 10 ksec in MCG-6-30-15 on the basis of an ASCA 
study. But this result has been recently questioned by Lee 
et al. (1999) and Chiang et al. (2000) who have studied 
with RXTE the spectral variabUity of MCG-6-30-15 and 
of NGC 5548, and have shown that the Fe Kq line flux 
is constant over time scales of 50 to 500 ksec, while the 
underlying continuum displays large flux and spectral vari- 
ations, and the proflle of the line is variable! The same re- 
sult seems to be observed in NGC 3615. Reynolds (2000) 
extended the work of Lee et al. (1999) on MCG-6-30-15, 
excluding a time delay between the line and the contin- 
uum in the range 0.5 to 50 ksec, and suggesting that the 



line flux remains constant on timcscalcs between 0.5 and 
500 ksec. From an analysis of the ASCA archive data of a 
sample of Seyfert nuclei. Weaver et al. (2001) found that 
in most cases changes in the line do not appear to track 
changes in the continuum. 

These problems indicate, that either the geometry of 
the accretion flow envisioned as an X-ray lamp above the 
flat disk plane is not correct, or the response of the mate- 
rial is not as expected from simple estimates, or both. 

In the present paper we concentrate on the second 
problem, i.e. on an accurate description of the X-ray re- 
processing by an optically thick material. Since we do not 
want to restrict ourselves to a particular model of accre- 
tion flow we discuss two basic general cases: we consider 
a medium at constant pressure (with the density gradient 
determined by the temperature profile and the radiation 
plus gas pressure gradient) and we compare this situation 
to the case of a medium at constant density. We assume a 
plane parallel geometry. In the next paper we will explore 
the infiuence of the geometry, with particular attention 
given to a quasi-spherical case. 

The radiative transfer computations arc done with the 
coupled codes titan and noar described in detail by Du- 
mont et al. (2000; hereafter DAG). We present a sequence 
of models characterized by a broad range of values of the 
ionization parameter defined at the surface. We also study 
the effects of other parameters involved, like the shape of 
the incident radiation flux. 

We show the broad band spectra, characterize their 
overall properties by several parameters frequently mea- 
sured observationally, like the index and the amount 
of reflection as seen in XSPEC analysis, and flnally we 
discuss the results comparing these quantities to typical 
parameters in various types of AGN. 

2. Radiation transfer 

The code titan was developed for the purpose of radia- 
tive transfer computations in an optically thick irradiated 
medium. Basically similar codes were developed by Ross 
& Fabian (1993), Nayakshin et al. (2000; hereafter NKK), 
Madej & Rozanska (2000). Other codes either do not cover 
the optical/UV energy range (Zycki ct al. 1994), or do not 
include hard X-ray irradiation (Hubeny et al. 2001), or 
they include hydrogen and helium only (Sincell & Krolik 
1997). 

The code titan, the code of Ross & Fabian (1993) 
and the NKK code based on xstar include the non-LTE 
effects while Madej & Rozanska (2000) still work in the 
LTE approximation. Both Ross & Fabian (1993) and NKK 
include the effect of Comptonization on the spectral dis- 
tribution of photons while titan considers only the ef- 
fect of Compton heating/cooling in energy balance. In 
order to obtain the Comptonized spectrum, the transfer 
code titan has to be supplemented with the Monte Carlo 
code NOAR of DAC. The code of Ross & Fabian (1993), 
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adopted by Ballantync, Ross & Fabian (2001), contains 
only a small number of basis atomic lines. The major dif- 
ference of TITAN with other codes (except that of Hubeny) 
is that it solves the line transfer instead of using the escape 
probability approximation. 

The detailed description of the code titan and exem- 
plary results for the temperature profile, ionization stage 
profile and the spectra in case of constant density of the 
medium are given by DAC and in Dumont & Collin (2001). 
The code was already used to reproduce the mean quasar 
spectra within the frame of the cloud model (Czerny & 
Dumont 1998), to model the iron line (after supplement- 
ing code results with NOAR: Abrassart 2000, Abrassart & 
Dumont 1998, 2000, Karas et al. 2000) and to model the 
reflection by an accretion disk (after coupling with the 
computations of the hydrostatic equilibrium: Rozanska et 
al. 2002). 

3. Results 

Since the dynamics of the formation of the hot and cold 
medium close to a black hole is not known it is necessary 
to reduce the problem to the consideration of a few basic 
possibilities for the description of the irradiated medium. 

If the cold medium consists of clumps of cold material 
formed through a violent disruption of the disk accretion 
flow, the clouds may not reach an equilibrium and may 
be represented as a constant density medium. However, if 
the clouds are illuminated and can survive long enough, 
the medium will approach pressure equilibrium and an as- 
sumption of constant pressure may be a better representa- 
tion of the medium, with the density gradient maintained 
due to the temperature and radiation plus gas pressure 
gradient. Finally, if the disk is not disrupted and it ex- 
tends down to the marginally stable orbit (and possibly 
below it) the irradiated upper layers of the disk will be 
mostly in hydrostatic equilibrium. The last possibility is 
characterized by the steepest density gradient. 

In the present paper we concentrate on the interme- 
diate case, i.e. the medium in pressure equilibrium. We 
assume a plane-parallel geometry, and the medium is ir- 
radiated by an isotropic incident flux. The irradiated lay- 
ers expand, and in the equilibrium situation a constant 
pressure is maintained throughout the entire zone. This 
pressure includes both the gas pressure and the radiation 
pressure. 

3.1. Constant pressure models 

We show flrst a set of models parameterized by the value of 

the ionization parameter ^ at the surface of the irradiated 
medium (z = 0). This parameter is deflned as 

£= ^ (1) 

where Lboi is the bolometric luminosity of the source lo- 
cated at a distance R and n{z = 0) is the number density 



at the surface of the medium. In our computations the 
density was fixed at the value n(0) = lO^^cm"^. 
Radiation pressure is locally calculated from 
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where is the mean intensity at a given depth computed 
at each step by the code. 

The value of the total pressure (i.e. radiation -|- 
gas pressure) is calculated when the temperature of the 
plasma has been determined by the thermal balance equa- 
tion, and this value is kept fixed within the medium thus 
determining the density profile. The model is iterated un- 
til convergence of the flux (and therefore of the radiation 
pressure) as in DAC. 

Our basic sequence of models is specified as follows. 

The shape of the incident radiation was assumed to be 
a power law, with the energy index a = 1.0, extending 
from 0.1 eV to 100 keV. 

The total hydrogen column in these computations is 
fixed at 3 X lO^^cm'^ (rrh » 24). 

The set of models with ^ = 300, 10^ 10^, 10^ and 10^ 
illustrates the effect of increasing the incident fiux from 
2.7 X 10^^ erg s~^ cm~^ to 8 x 10^^ erg s~^ cm~^ corre- 
spondingly. 

In all cases the irradiation results in the formation of 
a hot surface zone with an almost constant temperature. 
Initially the temperature decreases inwards rather slowly, 
with a rapid drop by 2 or 3 orders of magnitude below a 
certain depth. 

This rapid transition refiects the ionization instabil- 
ity which develops in an irradiated medium (Krolik et 
al. 1981), and it leads to a temperature discontinuity 
if the medium is supposed to be in hydrostatic equilib- 
rium (Begelman et al. 1983, Rozanska & Czerny 1996, 
Ko & Kallman 1994, NKK). Unique, continuous solu- 
tions are obtained only if energy transfer by conduction 
is included, in addition to radiative heating and cooling 
(Begelman & McKee 1990, Rozanska 1999). However, at 
present no computations include both the full radiative 
transfer and the conduction term; either radiative trans- 
fer is replaced by a parametric description of heating and 
cooling (Rozanska 1999) or the radiative transfer is com- 
puted carefully but the conduction is neglected (NKK, 
Ballantyne et al. 2001, Rozanska et al. 2002 and other 
papers presenting spectra). In the second approach the 
instability is dealt with, either by adopting a discontinu- 
ity at an arbitrary position somewhere inside the multiple 
solution zone, or the smooth solution is foimd at the ex- 
pense of a small pressure inversion in the multiple solution 
zone (see detailed discussion of this issue by Rozanska et 
al. 2002) . The conditions for the existence of multiple solu- 
tions under various approaches to line heating and cooling 
will be discussed by Coupe et al. (in preparation). In the 
present paper we adopt the continuous solution approach, 
as in Rozanska et al. (2002). 
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Fig. 1. The temperature (panel a), density (panel b) and pressure (panel c) profiles as a function of the optical depth 
in Thomson units, in a constant pressure medium, for an ionization parameter at the surface ^ — 300, 10'^, 10^, 10*^ and 
10^ (continuous, dotted, short-dashed, long-dashed and dash-dot line, correspondingly); the density at the surface is 
n(0) — lO^^cm"^, the column density is A^h = 3 x lO^^cm"^. The temperature profile as a function of the distance 
from the surface is shown in panel d. 



We see in Fig. |^ (panel a) that both the surface tem- 
perature and the optical thickness of the hot layer de- 
pend strongly on the incident flux. Only the outermost 
layers of t^s ~ 0.3, and only at high irradiation flux, reach 
the Inverse Compton temperature. Tic, corresponding to 
the incident spectrum (1.8 x 10^ K), mostly because of 
the contribution coming from the thermalized back side 
radiation. For ^ < 10'^ the surface temperature is much 
lower than Tic since the ionization of the upper layers is 
not high enough to make the Compton heating/cooling to 



dominate. The optical depth of the hot zone ranges from 
Th ~ 0.6 in the low flux case to Th ~ 2. Still lower ioniza- 
tion flux would lead to still lower optical depth of the hot 
skin and even lower surface temperature value. We do not 
show an example since the temperature of the back side 
of the slab approaches low values (T < 10* K) for which 
the code is not adapted. 

As we assume a constant pressure the density profile 
results from the computations (see Fig. 0, panel b). The 
density contrast between the illuminated and back sides 
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of the slab is very high since it reflects the temperature 
gradient as well as the considerable contribution from the 
radiation pressure which, in the case of high ^, dominates 
the gas pressure by several orders of magnitudes in the 
hot layers (see Fig. |l|, panel c). 

This density gradient means that, while the cold re- 
gion dominates the medium in terms of optical depth, the 
opposite is true if the spatial extension of the medium is 
considered. For example, for ^ ~ 10'^ the hot region is 
5 X 10^^ cm deep and the cold region is 10^'^ cm, while 
for ^ = 10^ the hot region is 8 x 10^^ cm and the cold 
region only 4 x 10^ cm! Of course the cold medium can be 
geometrically and optically thicker but we do not consider 
larger values of the column density since that would not 
change the reflected spectra. Therefore, apart from the 
temperature plots as function of the optical depth we also 
plot the temperature against the true distance from the 
irradiated surface (see Fig. |l|, panel d). 

The large relative geometrical size of the heated layer 
has an important consequence in the case of the special 
geometry of a quasi-spherical distribution of clouds irradi- 
ated by the central source, as discussed by CoUin-Souffrin 
et al. (1996). If the clouds have time to reach pressure equi- 
librium their irradiated parts would expand enormously 
and the medium will actually look more like clouds embed- 
ded in a warm bath, particularly for high irradiation flux. 
This will be particularly important when we will discuss 
the spherical geometry (Dumont et al., in preparation). 

As we mentioned earlier, the ionization of iron is com- 
plete only close to the surface, and only in the case of the 
highest irradiation flux. Even for ^ = 10® a large part of 
the hot layer is only partially ionized and populated with 
hydrogen and helium- like iron (see Fig. ||). The solution 
for ^ = 10^ is still mostly populated by FeXXV. Only in 
the case of the lowest value of the ionization parameter the 
contribution from highly ionized iron becomes negligible 
and the line is dominated by the 6.4 keV component. In 
all other cases the iron line consists always of a number of 
components at various energies which reflect the thermal 
stratification of the medium. 

The reflection spectra are given in Fig. ^. The ther- 
malized fraction of the incident radiation decreases with 
increasing ^ because the reflected fraction increases, due 
to the increased value of the temperature and of the opti- 
cal depth of the hot layer, and the temperature of the UV 
bump increases because the temperature of deep layers in- 
creases. The Big Blue Bump due to the thermalized radi- 
ation is rather narrow and it never extends to the observed 
soft X-ray band. The same effect occurs in the case of a 
medium in hydrostatic equilibrium (NKK, Rozahska et al. 
2002). In the far UV band the Lyman edge can be seen 
in emission (for low values of ^) or weakly in absorption 
(for ^ > 10^). This is simply due to the fact that the un- 
derlying continuum becomes more intense in this band as 
^ increases. Soft X-ray lines superimposed on a reflected 
continuum are clearly visible in all cases. The continuum 
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Fig. 2. The population of iron ions as a function of optical 
depth in a constant pressure medium, for ^ — 10'^ (upper 



panel) and 10^ 
3 x 1025cm-2. 



(lower panel); n(0) = lO^^cm , = 



itself clearly shows the signatures of the ionized reflection 
since it is relatively flat on vFi, plot and slightly rising 
below 2 keV, although not as strongly as in the constant 



density solutions (cf. Sect. 3.2). 

The iron Kq, line in all models consists of several 
components. The neutral component dominates only for 
^ = 300. In all other cases the line comes mostly from 
ionized iron ions. The helium-like iron contribution dom- 
inates for ^ = 10'^, and the role of the hydrogen-like iron 
increases with further increase of ^. This is not surpris- 
ing if we consider the stratification of iron ion populations 
(see Fig. H). The details of the spectra in this energy band 
determined from NOAR are shown in Fig. ^. Equivalent 
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Fig. 3. The reflection component from a constant pressure medium for ^ 
10"cm"3,iVH = 3 X lO^^cm"^^ Spectral resolution: 30. 



300, 10^ 10", 10^ and 10^; n(0) 



widths of the line components determined with respect to 
the total (i.e. reflected plus incident) radiation flux are 
given m Table 0. 

The line is broadened by Comptonization and a red 
tail of the line develops for large values of the ionization 
parameter. A narrow core, or rather several narrow com- 
ponents remain, and the effect of Comptonization can be 
seen as an additional broad underlying component. 

3.2. Comparison with a constant density case 

The reflection from the constant density medium was al- 
ready broadly discussed in the literature (Ross & Fabian 
1993, Zycki et al. 1994, Dumont & Collin 2001), including 



the Fe line properties (Matt et al. 1993, Zycki & Cz- 
erny 1994). However, the comparison is much easier if the 
computations are made with the same code and for the 
same sets of corresponding parameters. 

Therefore we show two examples, computed for an in- 
cident power law spectrum, the density equal to 10^^ cm~'^ 
everywhere in the medium and the total column density 
to 3 X 10^^ cm~^ as before. 

The temperature profile (compare Fig. |l| and |5|) is 
much more shallow that in the constant pressure case, 
the heated zone is effectively much broader. The overall 
spectrum (see Fig. |^) varies much more strongly with the 
incident flux than in the constant pressure case, and dis- 
plays a much broader "Big Blue Bump" , due to the pres- 
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Fig. 4. The total spectrum (reflection plus primary) in 
arbitrary units in the region of the iron line, computed 
with the NOAR code for a constant pressure medium; ^ ~ 
300,103,10^,105 and 10^ n(0) = 10"cm-3^7VH = 3 x 
lO^^cm-^. Resolution: 100. 



ence of layers with intermediate temperature. This strong, 
almost a power law, tail extends up to 2 keV thus forming 
a strong soft X-ray excess. 

Despite the more shallow temperature profile the iron 
line in the constant density case is usually dominated 
by a single component. This is due to the fact that the 
zone where the temperature traces the ionization state is 
Thomson thick and the line is effectively dominated by 
the ionization state of the outer layers thus giving a cold 
line for low ^, and 6.7 - 6.9 keV lines for higher irradia- 
tion flux. Any particular ionization stage of the iron can 
be modeled by adjusting the incident flux, which is not 
the case either in constant pressure computations or in 
solutions in hydrostatic equilibrium, as noticed by NKK. 

In the case of high ionization parameter — 10^) 
the line again is broadened by Compton downscattering 
(see Fig. ^ and the effect is stronger in the constant den- 
sity medium than in the constant pressure case (see also 
Abrassart & Dumont 2000 and Abrassart 2000). In the 
case of low ^ the line is dominated by the narrow cold 
component but we also see a strong absorption edge. Such 
a feature is only barely seen from the constant pressure 
medium. 
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Fig. 6. The reflection component from a constant density 
medium, for ^ = 300, iVn = SxlO^^cm^^ ^ ^nd for ^ = 10^; 
iVn = lO^^cm"^; n(0) = lO^cm'^. Resolution: 30. 



3.3. Comparison with hydrostatic equilibrium 

Two examples of reprocessing by an irradiated accretion 
disks were presented by Rozahska et al. (2002). The pa- 
rameterization of such models is through the accretion 
rate, and the pressure increases towards the disk equa- 
torial plane. However, we can compare two models with 
similar total pressure at the surface. 
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Fig. 7. The total spectrum (reflection plus primary) in ar- 
bitrary units in the region of the iron line, computed with 
the NOAR code for a constant density medium; ^ — 300, 
iVn = 3 X lO^^cm-^ , and for ^ = 10^ TVh = 
n(0) = 10"cm-3. Resolution: 100. 



In the model computed for a black hole mass M = 
IO^Mq, an accretion rate m = 0.03AfEdd and a distance 
10i?schw, the gas pressure is of the order of 2 x 10^ g cm~^ 
s~^ in the surface layers, so it approximately corresponds 
to our case ^ — 10^ of constant pressure models. 

The surface temperature is slightly lower in the con- 
stant pressure model since the slope of the incident radi- 
ation is different in the two models (a — 0.9 was adopted 
in the hydrostatic equilibrium model, while a = 1.0 in the 
present computations). The overall shape of the reflected 
spectrum in X-ray band is similar, apart from the effect 
caused by the difference in the adopted incident flux in the 
two computations. The iron Kq, line in both cases shows 
all three components (6.4 keV, 6.7 keV and 6.9 keV), but 
the neutral component is relatively stronger in the case of 
hydrostatic equilibrium. The Big Blue Bump peak is also 
relatively stronger in the case of hydrostatic equilibrium. 
This is partially due to an additional energy generation in 
the accretion disk itself and partially due to a much thicker 
hot skin in the case of constant pressure medium (~ 0.7 
instead of ~ 0.1) which increases the reflected fraction at 
the expense of the reprocessed one. Increased Thomson 
thickness of the radiatively heated zone accounts also for 
the difference in the value of the equivalent width of the 
iron Kq, line - the line from the constant pressure model 



is by 50% weaker than in the hydrostatic equilibrium case 
(66 eV vs. 92 eV). 

Similar differences are seen between the irradiated disk 
model rh = 0.3 (hydrostatic equilibrium) and the (con- 
stant pressure) model with ^ = 10^, both having a to- 
tal surface pressure of order of 2 x 10^. Also in this case 
the skin is thicker in the constant pressure case (rxh = 1 
vs. TTh = 0.3), and consequently the equivalent width of 
the Fe Kq, line is lower (80 eV in hydrostatic equilibrium 
case and 53 eV in constant pressure case), and the domi- 
nance of the 6.7 keV component in the line proflle is much 
stronger in the constant pressure solution. 

3.4- The effect of the spectral shape of incident radiation 

3.4.1. Power law slope 

The slope of the incident radiation is important not only 
because of its contribution to the observed spectrum. It 
also modifles the spectral shape of the reflected component 
by affecting the value of the Compton temperature and the 
thickness of the hot layer. 

In constant density computations the incident slope 
determines the reflected slope in the soft part of the spec- 
trum (see Fig. 6 of Czerny & Dumont 1998) - a steeper 
incident flux results in a proportionally steeper reflected 
component. The properties of the iron line do not change 
dramatically. 

A medium in hydrostatic equilibrium reacts differently 
to such changes (NKK): a harder spectrum results in a 
completely ionized outer skin and a relatively faint but 
neutral iron line produced in the deeper un-ionized zone, 
while a softer spectrum gives only a partially ionized skin 
and a strong iron line coming from ionized iron ions in 
that region. 

The reaction of a medium at constant pressure is in- 
termediate between these two cases (see Fig. ^ . When the 
radiation spectrum is very hard (a — 0.7) the dominant 
contribution to the line comes indeed from the deep lay- 
ers and the 6.4 keV component dominates. Intermediate 
slopes enhance the contribution from helium-like iron and 
both line components (6.4 and 6.7 keV) are clearly seen. 
However, when the spectrum is very soft (a = 1.3), a cold 
iron line dominates again since the gas is not heated to a 
high enough temperature (even at the surface) for highly 
ionized iron to contribute. 

3.4.2. The spectrum of Laor et al. (1997) 

A hot (thermal) plasma is a source of a power law ra- 
diation which is formed due to Compton upscattering of 
the soft photons by energetic electrons. A cold layer ther- 
malizes a fraction of the radiation but nevertheless, in the 
plane-parallel geometry, the X-ray primary emission dom- 
inates the total spectrum. This is not a good representa- 
tion of high accretion rate objects like quasars or NLSl 
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Fig. 8. The dependence of the reflection component from a constant pressure medium on the slope of the incident 
radiation spectrum a equal 0.7, 0.9, 1.0 and 1.3, for ^ = 10'^; n(0) = lO^^cm"'', A'^h = 3 x lO^^cm"^. Spectra are shifted 
for better illustration. Resolution: 30. 



galaxies (where the soft thermal component dominates), 
which means that either the geometry is different, or most 
of the energy is dissipated within the cold phase at larger 
distances from the black hole. In this latter case the con- 
sidered innermost cold layer will see not only the hard 
X-ray source but it will be also affected by the mean ra- 
diation field. 

To illustrate this situation we consider the case of a 
layer irradiated by a source with a spectral distribution 
given by the mean spectral shape of quasars, as deter- 
mined by Laor et al. (1997) from a large quasar sam- 
ple. This radiation field is much softer than our assumed 
power law, since it is strongly dominated by the Big Blue 
Bump, and the corresponding Inverse Compton tempera- 



ture is lower than for a power law spectrum with a = 1.0 
(2.4 X 10'' K and 1.9 x 10'' K, correspondingly). 

The reflected spectrum in the case of a constant pres- 
sure medium is shown in Fig. ^. It traces to some extent 
the incident spectrum, giving a strong soft X-ray excess, 
unlike in the case of an incident power law spectrum. The 
overall shape is much more similar to the case of con- 
stant density computations, with a power law spectrum. 
A strong soft X-ray excess is clearly seen below 2 keV. 

In the Kq, line profile a strong contribution from 6.4 
keV component is always seen, although for ^ < 10^, the 
6.7 keV component is slightly stronger (see Fig. |l^. A 
model with a power law incident radiation and a compa- 
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Fig. 9. The reflection component from a constant pressure 
medium for the incident radiation spectrum given by Laor 
et al. (1997), for ^ = 10^ and lO'^. Parameters: n(0) = 
10"cm-3,iVH = 3 X 1025cm-2. Resolution: 30. 
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rable hard X-ray flux does not show such a strong neutral 
component. 

3.5. Total column density 

Our computations were dehberately made for a Thomson 
thick medium. However, an interesting question is, how 
strong the effect of the column density is for a moderately 
thick medium, and what optical thickness of the medium 
is needed to thermalize the incident radiation. 

In the case of a constant density medium the photons 
penetrate quite deeply into the medium if the irradiation is 
strong. We saw in Fig. ^ that for ^ = 10^ the temperature 
of the back side was still quite high even for the very large 
column density adopted in this computation (10^^ cm~^). 

However, in constant pressure medium the thermaliza- 
tion is much more efficient for the same Thomson optical 
depth due to the fact that the density is much higher in the 
cold region, and therefore coUisional cooling is much more 
efficient. We give two examples of the reflection spectra 
calculated for a moderate value of the ionization parame- 
ter = 10'^) and two values of the column density (cor- 
responding to TTh = 20 and rxh = 2; Fig. |l^). Even if the 
total optical thickness of the zone is equal to 2 the tem- 
perature at the back side of the layer is low. The resulting 
spectra are almost indistinguishable. 

3.6. Density 

Parameterization of the incident flux with the ionization 
parameter ^ has the advantage of reducing significantly 



Fig. 10. The total spectrum (reflected plus primary) in 
the region of the iron line for two sets of models, with 
similar normalization of the incident flux at 2 keV but two 
shapes of the incident radiation spectra. First set: power 
law spectrum with ^ = 300 and the Laor et al. (1997) 
model with ^ = 10'^. Second set: power law spectrum with 
^ = lO'* and the Laor et al. (1997) spectrum with ^ = 
10^ from a constant pressure medium. Other parameters: 



n(0) = 10 



3 X lO^^cm-^. Resolution: 100. 



the dependence on the medium density. In constant den- 
sity computations this approach was particularly succes- 
fuU and the reflection spectra do not depend significantly 
on the density in a broad range of parameters (from 
n — 10^°cm~'^ up to n = 10"'^'*cm"'^). Some differences are 
seen in the low energy band (IR and optical) due to the 
influence of free-free emissivity and the reflection hump 
above w 3 keV is more prominent in the low density case 

(Fig. li). 

Radiation transfer computations performed in con- 
stant pressure medium are much more sensitive to the 
density itself due to the density gradient. When the sur- 
face density is changed the radiation to gas pressure ratio 
changes at the surface and the new density profile as a 
function of the optical depth changes, and it is not just 
a rescaled version of the old profile by a constant value. 
This induces some differences in the broad band spectrum 



as well as in details of the spectral features (see Fig. 13). 
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Fig. 11. The reflection component from a constant pres- 
sure medium for ^ — 10"^; n(0) — 10^^ cm~^, and two 
values of the column density: A'h = 3 x lO^^cm"^ and 
A^H = 3 X lO^'^cm"^. The spectra are shifted for a better 
illustration. Resolution: 30. 
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Fig. 13. The reflection component from a constant pres- 
sure medium for the power law incident radiation spec- 
trum for ^ — 10"^ and two values of the number den- 
sity; n(0) = 10"cm-3 and n(0) = IQi^cm-^; iVjj = 
3 X 1025cm-2. 




Fig. 12. The reflection component from a constant density 
medium for the power law incident radiation spectrum for 
^ = 10^ and two values of the number density; n(0) = 
10i°cm-3 and n(0) = IQi^'cm-^; A^h = 3 x lO^^cm-^. The 
spectra are shifted for a better illustration. Resolution: 30. 



3. 7. Summary of models 

In order to facilitate the comparison between the vari- 
ous models presented here and the observed properties 
of the AGN spectra we determine for those models sev- 
eral broadly used quantities. The results are presented in 
Table |l| and they include more models than just those 
presented in the previous subsections. 

Since the computations are made in plane-parallel ge- 
ometry the final spectrum consists of the reflected spec- 



trum plus the incident spectrum, without any additional 
weight. 

We introduce the aox parameter as the mean slope cal- 
culated between 2500 A and 2 keV. We give the widths of 
the components of the Fe Kq, line, measured with respect 
to the observed continuum (i.e. reflected plus incident) 
which result directly from our computations. 

Since the observed spectra are usually analyzed using 
the XSPEC software (and the standard spectral compo- 
nents provided with it) we follow other authors presenting 
their theoretical results (Done & Nayakshin 2001, Ballan- 
tyne et al. 2001) and we analyze our spectra as we would 
analyze an observed spectrum. 

We create artificial data using the Jake command in 
XSPEC for ASCA satellite and assuming the count rate 
at the typical level expected from bright AGN sources with 
a long exposure. 

We fit such data using pexrav model (Magdziarz & 
Zdziarski 1995) to represent the incident power law and 
the reflected continuum, and we add GAUSSIAN to account 
for the iron line. These fits provide us again with the 
normalization of the reflected component, i?xsPEC, and 
the photon index of the incident radiation, T. In princi- 
ple i?xsPEC should be equal to 1, while F should be equal 
to the photon index adopted in our model computations. 
However, the shape of the reflected component assumed in 
PEXRAV is rather different from the one from our compu- 
tations of radiative transfer and so both i?xsPEC and T are 
affected. The results for various models are summarized in 
Table |l|. We see that low ^ models show a reflection com- 
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ponent larger than 1, while reflection from highly ionized 
material give values significantly lower than 1. This effect 
is much stronger in the constant density models which 
are characterized by a reflection coefficient of the order of 
unity for a rather narrow range of values of the ioniza- 
tion parameter. Fits to models with the Laor et al. (1997) 
incident spectrum gave formally very large values of the 
reflection coefficient which was caused by the presence of 
the strong soft X-ray excess in those spectra models. We 
did not add any additional component during the fitting 
to account for this excess which leads to an artificially 
steep 'incident radiation' and a large reflection. 

We also considered the XSPEC pexriv model 
(Magdziarz & Zdziarski 1995) of ionized reflection. How- 
ever, since the shape of the reflection component is differ- 
ent in our computations and in PEXRiv, this model did 
not provide a better description of our model and the 
returned value of the ionization parameter was usually 
rather low independently from the true value characteriz- 
ing the model. Similar problems were found by Ballantyne 
et al. (2001). 



4. Discussion 

X-ray reprocessing in plane-parallel geometry under the 
assumption of the constant density has been studied in a 
number of papers (Ross & Fabian 1993, Zycki et al. 1994, 
DAC; see Dumont & Collin 2001). Similar computations 
wore also performed for a stratified medium in hydrostatic 
equilibrium (Raymond 1993, NKK, Ballantyne et al. 2001, 
Rozanska et al. 2002). 

In the present paper we consider an intermediate case: 

the medium under constant pressure. It may represent a 
strongly flattened distribution of clumpy accretion flow. 

Our model, as most of the constant density or hydro- 
static models, is clearly oversimplified since we assume 
that the entire flow is characterized by a single value of 
the ionization parameter, which may not be the case in 
a real accretion flow. We also did not include here the 
dissipation which may take place in the cold phase of an 
accretion flow and consequently enhance the contribution 
of the Big Blue Bump to the overall spectrum. This last 
effect is only simulated to some extent by considering the 
shape of the incident radiation flux in the form of the mean 
quasar spectrum as determined by Laor et al. (1997). 

Nevertheless, the results illustrate certain important 
trends which should help to analyse the conditions in the 
emitting plasma. 

Below we compare the results obtained from the con- 
stant pressure medium with the observed trends in the 
broad band spectra of AGN. 



4.1. Observed properties in quasar and Seyfert 

galaxies samples 

The observed broad band spectral index aox measured be- 
tween 2500 A and 2 keV was shown for a sample of Seyfert 
galaxies by Walter & Fink (1993) and for quasars by Green 
(1996). Both samples showed a large dispersion of this 
quantity, which could be partly explained by reddening 
(e.g. Wilkes et al. 1999, Laor & Brandt 2002). The aver- 
age values are about 1.3 for Seyferts and 1.5 for quasars, 
with most Seyferts having values between 1 and 1.6, and 
most quasars having values between 1.2 and 1.8. 

Our constant pressure models with an incident power 
law radiation slope equal to 1 have on average somewhat 
lower values of the aox, around 1.2, with typical values 
between 1 and 1.3. Since the observed power law index in 
Seyfert galaxies is on average not steeper than the value 
adopted in our computations it shows that either there is 
some dissipation within the cold medium which slightly 
contributes to the Big Blue Bump, or the plane-parallel 
geometry adopted in the computations is not correct (for 
spherical accretion models see CoUin-Souffrin et al. 1996, 
Czerny & Dumont 1998, and Malzac 2001). The effect is 
even stronger in the case of quasars. Not surprisingly, the 
models with an incident flux based on averaged quasar 
spectra reproduce the observed values of aox for those 
models. 

Hydrostatic equilibrium models give for aox re- 
sults comparable to the constant pressure models (see 
Rozanska et al. 2002), with values slightly lower than seen 
in Seyfert galaxies, and lower than in quasars. Constant 
density models give still lower values (~ 1) and a bigger 
discrepancy between the simple prediction and the data. 

4-2. Amplitude of the reflection component 

The observed range of the reflection amplitudes for Seyfert 
galaxies (~ 0.3 — 1.5; Zdziarski et al. 1999) is rougly con- 
sistent with the values obtained from constant pressure 
models (0.6 - 1.8). Models with constant density easily 
produce values of the amplitude of reflection either much 
larger or much lower than observed. 

However, such a simple model does not explain the ob- 
served correlation between the slope of the incident radia- 
tion and the amoimt of reflection. The observational data 
show a strong increase of the amount of reflection with 
the steepening of the incident radiation. Our set of mod- 
els shows an opposite trend: a model with a steep (soft) 
incident radiation flux (a = 1.3) shows a fitted amount of 
reflection much lower than a model with a flat (hard) inci- 
dent radiation (a = 0.7). However, our computations were 
done for a single value of the ionization parameter ^ = 10^ 
which might affect the results. We could not repeat the 
same computations for a practically neutral medium since 
for the moment our codes do not allow to solve the cou- 
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Table 1. Summary of model properties. 



No. 


C incid*" 






-'-'max 


'-''OX 


(6.4) 


(6.7) 


(6.9) 


ryh 

-KXSPEC 




constant pressure 


1 


3 X 10^ 


1.0 


1 X 10" 


3 X 10^^ 


5.21 


1.31 


18.0 


3.4 


0.0 


1.53 ±0.17 


2.09 ± 0.05 


2 


1 X 10^ 


1.0 


1 X 10" 


3 X 10^^ 


6.45 


1.25 


14.9 


96.7 


5.7 


1.82 ±0.10 


2.11 ±0.03 


3 


1 X 10"* 


1.0 


1 X 10" 


3 X 10^^ 


8.39 


1.13 


4.5 


42.8 


19.2 


1.14 ±0.05 


2.06 ±0.02 


4 


1 X 10^ 


1.0 


1 X 10" 


3 X lO^"' 


12.62 


1.04 


5.8 


37.2 


9.7 


0.88 ±0.08 


2.05 ±0.02 


5 


1 X 10'^ 


1.0 


1 X 10" 


3 X 10^^ 


22.68 


1.01 


2.9 


18.8 


8.8 


0.58 ± 0.02 


2.05 ± 0.02 


6 


1 X 10^ 


1.0 


1 X 10" 


3 X lO^-* 


6.39 


1.24 


74.9 


102.6 


4.0 


1.61 ± 0.09 


2.10 ± 0.03 


7 


1 X 10^ 


0.7 


1 X 10" 


3 X 10^^ 


2.37 


1.16 


16.0 


78.7 


4.2 


1.71 ±0.07 


1.81 ±0.02 


8 


1 X 10^ 


1.3 


1 X 10" 


3 X 10^^ 


5.97 


1.58 


33.9 


0.0 


0.0 


0.87 ±0.04 


2.06 ±0.03 


9 


1 X 10^ 


1.0 


1 X 10^^ 


3 X 10^^ 


19.64 


1.11 


87.0 


99.8 


2.1 


1.83 ±0.03 


2.11 ±0.02 


10 


1 X 10^ 


Laor 


1 X 10" 


3 X 10^^ 


11.30 


1.54 


70.0 


91.1 


13.0 


17.21 ± 0.45 


2.77 ±0.05 


11 


1 X 10^ 


Laor 


1 X 10" 


3 X 10^^ 


6.50 


1.66 


114.5 


0.0 


0.0 


17.95 ± 0.55 


2.76 ± 0.05 


constant density 


12 


3 X 10^ 


1.0 


1 X 10" 


3 X 10^" 


4.97 


1.15 


19.9 


4.0 


0.0 


2.19 ±0.18 


2.12 ±0.01 


13 


1 X 10^ 


1.0 


1 X 10" 


1 X 10^® 


368.6 


0.99 


0.0 


0.0 


1.0 


0.15 ±0.06 


2.02 ± 0.01 


14 


1 X 10'"* 


1.0 


1 X 10^^ 


3 X 10^"^ 


10.4 


1.18 


0.0 


105.9 


4.0 


1.32 ±0.08 


2.03 ±0.03 


15 


3 X 10^ 


Laor 


1 X lO^-* 


3 X 10^^ 


9.86 


1.55 


89.8 


0.0 


0.0 


11.13 ±0.01 


2.78 ± 0.02 


16 


2 X 10® 


Laor 


1 X 10" 


3 X lO^-* 


11.11 


1.50 


0.0 


0.0 


7.1 


14.34 ± 0.01 


2.69 ± 0.02 



" - ionization parameter 

- incident spectrum shape (either AGN mean spectrum after Laor et al. (1997), or a power law with an energy index (a) 

- number density at the surface 

- total column density of the medium 

- energy at which the spectrum has its maximum at the EUVE range in [eV] 
^ - effective spectral index measured between 2500 A and 2 keV 

- amount of reflection as determined from pexrav XSPEC model 
^ - X-ray slope as determined from pexrav XSPEC model 



pled equations of radiative transfer and energy balance if 
the temperature is too low. 

Hydrostatic equilibrium models give lower values of 

the reflection component 0.5 for the incident spectrum 
slopca 1) than constant pressure medium, and they 
show a clear correlation between the slope of the inci- 
dent radiation and the amount of reflection, although the 
range of values is somewhat smaller than observed (Done 
& Nayakshin 2001). 

Constant density models exhibit the biggest range in 
the amount of reflection (from to 2), rather too large in 
comparison with the data. 

4-3. The Fe line diagnostics 

The structure and the broadening of the iron line seen 

in the data is still under discussion and the situation is 
not clear. The broad line profile seen in the ASCA data 
of the source MCG-6-15-30 was confirmed in subsequent 
observations. However, the shape of the iron line in a typ- 
ical Seyfert galaxy may be narrower than inferred initially 
(compare the results of Nandra et al. 1997 and Lubinski 
& Zdziarski 2001 based on ASCA data). Recent Chan- 
dra data point towards a very narrow neutral component 



(e.g. Yaqoob et al. 2001a, Sambruna et al. 2001, Kaspi et 
al. 2001) but it may be superimposed on a broader com- 
ponent originating closer to the central object (Yaqoob 
et al. 2001b). The line is therefore most probably multi- 
component but the components have different shapes and 
originate in different places. 

Our analysis show that the iron line originating from 
a single region is expected to be multicomponont if the 
medium is at constant pressure when the slope of the 
incident radiation spectrum is a 1, and the medium 
is not neutral (i.e. the ionization parameter larger than 
~ 200). The multicomponent structure originates in the 
same medium, so if kinematical broadening is essential, it 
should broaden all components similarly. 

A similar multicomponent line is present in the spectra 

of the hydrostatic models of Rozanska et al. (2002). Com- 
putations of NKK rather show a single neutral component 
for a ~ 1. 

A single component, either neutral or ionized, is in- 
stead characteristic for a constant density medium. 

There is some level of intrinsic broadening of the line 
profile due to Comptonization within the reprocessing 

medium, as suggested by Czerny et al. (1991) but the ef- 
fect is not very strong in the case of a constant pressure 
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medium or a medium in hydrostatic equilibrium. A con- 
stant density medium broaden the hne more effectively, 
due to the larger optical depth of the hot skin. 

At present, the observations are not accurate enough to 
differentiate between the constant pressure and the con- 
stant density medium, if the narrow neutral component 
comes from the outer parts of the accretion flow (BLR 
and/or dusty torus) and the broad component is barely 
seen and not well constrained. The study of variability 
could help to locate the emission regions of the different 
components, and subsequently help to find an answer to 
the question of hydrostatic equilibrium. 
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